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Abstract 

The purpose of this paper is the fundamental theory of the non-uniform motion 
of dislocations in two and three space-dimensions. We investigate the non-uniform 
motion of an arbitrary distribution of dislocations, a dislocation loop and straight 
dislocations in infinite media using the theory of incompatible elasto dynamics. The 
equations of motion are derived for non- uniformly moving dislocations. The re- 
tarded elastic fields produced by a distribution of dislocations and the retarded 
dislocation tensor potentials are determined. New fundamental key-formulae for 
the dynamics of dislocations are derived (Jefimenko type and Heaviside-Feynman 
type equations of dislocations). In addition, exact closed- form solutions of the 
elastic fields produced by a dislocation loop are calculated as retarded line inte- 
gral expressions for subsonic motion. The fields of the elastic velocity and elastic 
distortion surrounding the arbitrarily moving dislocation loop are given explicitly 
in terms of the so-called three-dimensional elastodynamic Lienard-Wiechert tensor 
potentials. The two-dimensional elastodynamic Lienard-Wiechert tensor potentials 
and the near-field approximation of the elastic fields for straight dislocations are 
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calculated. The singularities of the near-fields of accelerating screw and edge dislo- 
cations are determined. 

Keywords: dislocation dynamics; non-uniform motion; dislocation loop; elastody- 
namics; radiation; retarded fields; near-fields; singularities. 



1 Introduction 

The investigation of the non-uniform motion of dislocations is an important interdisci- 
plinary research field. It has attracted the attention of scientists from several different 
fields such as applied mathematics, solid state physics, material science, continuum me- 
chanics, and seismology (see, e.g, [1, 2, 3]). Typical problems of elastodynamics of dis- 
locations are the determination of elastic fields produced by the non-uniform motion of 
straight dislocations and dislocation loops. 

The theory of the non-uniform motion of Volterra dislocations has a long history 
starting with the famous and well-known papers of Eshelby [4, 5]. Using an electro- 
magnetic analogy, Eshelby [5] found the elastic fields of a non-uniformly moving screw 
dislocation. The problem of the non-uniform motion of a gliding edge dislocation was 
solved by Kiusalaas and Mura [6]. Kiusalaas and Mura [6] gave the solution in terms 
of 'stress (or potential) functions' for the velocity and elastic distortion fields (see also 
[7, 1]). Lazar [8] has written a systematic review paper about the non-uniform motion of 
straight dislocations including the solution of a non-uniformly climbing edge dislocation. 

The behaviour of a straight dislocation is somehow particular, because at any time 
the fields are determined not only by the instantaneous values, but also by the values 
in the past [4, 5]. As Eshelby [4] succinctly put it: 'The dislocation is haunted by its 
past'. For that reason, all the solutions of the elastic fields of non-uniformly moving 
straight dislocations are given in the form as time integrals and show an afterglow. Due 
to the afterglow, Huygens' principle is not valid in two dimensions (see, e.g., [9, 10, 11]). 
In general, a two-dimensional wave-motion possesses a 'tail'. The wave motion in two- 
dimensions and the afterglow effect are discussed more in detail by Baker and Copson 
[12], Barton [13], and L 

The non-uniform motion of dislocation loops was also studied. The solution of non- 
uniformly moving dislocation loops was first formulated by Mura [14] in terms of the 
three-dimensional elastodynamic Green tensor as double integrals over the loop curve 
and time (see also [15]). The mathematical formulation of moving dislocation loops and 
analogous double integral presentations were also given by Kossecka [16] and Kossecka 
and dcWit [17. 18]. It was pointed out by Markcnscoff [19] (see also [20, 21]) that some 
care is necessary in the calculation of the elastic fields produced by non-uniformly moving 
dislocations in order to avoid non-integrable singularities. Markenscoff [19] showed that 
the general expressions for the velocity and elastic distortion fields of dislocations given 
by Mura [14, 15] are not free of non-intcgrablc singularities. The reason is that the 
integration and differentiation cannot be changed in some cases. 

In general, the problem of moving dislocation loops is a three-dimensional problem in 
space. It is well-known that Huygens' principle is only valid for odd space dimensions: 



2 



> 3 [9, 10, 11]. For that reason, there is no afterglow in three dimensions. In standard 
books on electrodynamics (see, e.g., [22, 23]), the electric scalar and magnetic vector 
potentials of a non-uniformly moving point charge which are the famous Lienard-Wiechert 
potentials [24, 25] can be found. It is quite surprising that nothing has been done in this 
direction in the elastodynamics of moving dislocation loops up to now. In elastodynamics, 
only the retarded potentials were given for the waves produced by body forces, using 
the Helmholtz decomposition (see, e.g., [10, 26]). A more general expression for the 
retarded potential in elastodynamics was given by Hudson [27]. Though Emil Wiechert 
was a director of the geophysical laboratory in Gottingen, he formulated his theory for 
electrodynamics [25], but not for elastodynamics, with which he was certainly familiar. In 
this paper, we investigate some fundamental problems of the elastodynamical theory of 
dislocations in analogy to the electromagnetic field theory. Especially, we want to develop 
the Lienard-Wiechert potentials for the elastodynamics and to apply them to dislocation 
loops and straight dislocations. 

This paper is organized as follows. In Section 2, the framework of incompatible elasto- 
dynamics and the equations of motion of dislocations are presented. In Section 3, the equa- 
tions of motion for an arbitrary three-dimensional distribution of dislocations are solved, 
using the three-dimensional elastodynamic Green tensor. As a result the retarded elastic 
fields are given. The non-uniform motion of a closed dislocation loop is studied in Sec- 
tion 4. The so-called elastodynamic Lienard-Wiechert tensor potentials are determined. 
The elastic fields are given in terms of so-called Lienard-Wiechert tensor potentials. The 
static limit of the elastic fields of the non-uniformly moving dislocation loop is given in 
Section 5. In Section 6, the two-dimensional elastodynamic Lienard-Wiechert tensor po- 
tentials and the elastic fields of non-uniformly moving straight dislocations are calculated. 
In Section 7, the near-field approximation of accelerating screw and edge dislocations are 
calculated. In a straightforward manner, the 1 /-R-singularity and a logarithmic singular- 
ity associated with the acceleration of the dislocation are found. The relation between the 
elastodynamic Lienard-Wiechert tensor potentials and Mura's dislocation tensor poten- 
tials is presented in Section 8. The retarded dislocation tensor potentials and the proper 
elastodynamic Lienard-Wiechert tensor potentials of a dislocation loop are determined. 
In Section 9, the conclusions are given. 



In this section, we derive the equations of motion for dislocations in the framework of 
incompatible elastodynamics (see, e.g., [14, 15, 16, 28, 29, 8]). An unbounded, isotropic, 
homogeneous, hnearly elastic sohd is considered. In the elasticity theory of self-stresses 
the equilibrium condition is^ 



where p and a are the linear momentum vector and the force stress tensor, respectively. 
In the incompatible linear elasticity, the momentum vector p and the stress tensor <t can 



2 The equations of motion of dislocations 



Pi 1 




^We use the usual notation := dkl^ij and $ij := c^tftj- 
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be expressed in terms of the incompatible elastic velocity (particle velocity) vector v and 
the incompatible elastic distortion tensor /3 by means of the two constitutive relations 

Pi^pvi, (2) 

— Cijkl Pkl , (3) 

where p denotes the mass density and Cijki is the tensor of elastic moduh. The tensor 
Cijkl is characterized by the symmetry properties 

Cijkl = Cjiki = Cijik = Ckiij ■ (4) 

If the constitutive relations (2) and (3) are substituted in Eq. (1), the equilibrium condition 
expressed in terms of the elastic fields v and /3 is obtained 

pvi — CijkiPkij ■ (5) 

The presence of dislocations makes the elastic fields incompatible which means that they 
are not anymore simple gradients or time derivatives of a displacement vector u. Unlike 
the displacement field and the plastic fields, the elastic fields are physical state quantities 
of dislocations. This is one reason, why we deal only with the calculation of the elastic 
fields in this paper. 

Other important tensor fields of dislocations are the dislocation density and dislocation 
current tensors (e.g. [28, 30]). The dislocation density tensor T and the dislocation current 
tensor I are defined by (see also [31, 30]) 

Tijk PikJ Pij,k ) (6) 

The dislocation current tensor / was originally introduced by Kosevich [32, 33] and 
Hollander [34, 35] (see also [28, 29, 36, 37]). The dislocation current tensor (7) is the dif- 
ference of two pieces: the time derivative of the elastic distortion, and the elastic velocity 
gradient. It is noted that T^jk — —T^kj- Both T and I have nine independent compo- 
nents. A field theoretical justification for the structure of T and I was given by Lazar 
and Anastassiadis [38] (see also [30]). Moreover, they fulfill the Bianchi identities (see 
also [39, 40]) 

^jki Tijk,i = , (8) 

-^ij'fc ~l~ ^i3,k -^ik-ij i (9) 

which are 'conservation' laws. Here ejki denotes the Levi-Civita tensor. Eq. (8) states 
that dislocations cannot end inside the medium and Eq. (9) means that the time evolution 
of the dislocation density tensor T is determined by the 'curl' of the dislocation current 
tensor I. The tensors T and / may describe single straight dislocations, dislocation loops 
and an arbitrary distribution of dislocations. 

Alternatively, we may rewrite the tensor T as (dual) tensor of rank two 

Oiij = - ^jkiTiki — ejkiPii,k (10) 
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with the inverse relation 



Tiki — ^jklOiij . (11) 

The tensor a is the usual dislocation density tensor^. Then the Bianchi identities (con- 
servation laws) (8) and (9) simplify to (see also [28, 44, 29, 45, 36]) 

«yj = 0, (12) 
oiij + ^jkilik,i = . (13) 

Now we derive separated equations for the elastic fields /3 and v as equations of motion. 
If we differentiate the Eq. (5) with respect to and use Eqs. (6) and (7) to eliminate 
V, we get the equation of motion for the incompatible elastic distortion tensor (3 (see 
also [16, 36, 8]) 

P Pirn Cijkl Pkmjl Cijkl TkmlJ plim i (1^) 

where the dislocation density and the dislocation current tensors are the sources. Eq. (14) 
is a tensorial Navier equation for (3. Similarly, if we perform the differentiation of Eq. (5) 
with respect to time and use Eq. (7) to eliminate /3, we obtain the equation of motion for 
the incompatible elastic velocity vector v (see also [16, 36, 28, 8]) 

pvi — Cijkl Vk,ji — Cijkl hi,j , (15) 

where the dislocation current tensor is the source term. Eq. (15) is a vectorial Navier 
equation for v. We want to note that also Rogula [46] obtained equations of the form 
(14) and (15) for dislocations in a pseudo-continuum 



3 The retarded elastic fields 

In this section, we calculate the retarded elastic fields produced by a three-dimensional 
distribution of dislocations. The solutions of Eqs. (14) and (15) can be represented as 
convolution integrals [14, 16, 18, 8]. For an unbounded medium and under the assumption 
of zero initial conditions, which means that /3(r,to) and v(r,to) and their first time 
derivatives are zero for ^ — oo, the solutions of /3 and v can be represented as 

/t roc 
/ Cjkin Gij (r -r',t- t') Ti^nir', t') dr' dt' 
■oo J —oo 

/t POO 
/ p dj {r -r\t- t') I,m{r\ f) dr' dt' (16) 
-oo J —OO 

^In the literature, the notations of the dislocation density tensor and the dislocation current tensor are 
not unique: Q;(Lazar [30]) = Q;(Kossecka [16, 41]) = Q;''"(Kossecka and deWit [17]) = a'^ (Kroner [42]) = 
aT(deWit [43]) = -Q;(Teodosiu [36]) = -a'r(Kosevich [28]) and /(Lazar [30]) = /(Kossecka [16]) = 
-/"^(Kossecka and deWit [17]) = -/(Teodosiu [36]) = /"^(Kosevich [28]). 
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and 



/t POO 
/ Cjkim G^j{r -r\t- t') Iim{r', t') dr' dt' . (17) 
-oo J —oo 

It is obvious that dislocations act as a source of the elastodynamic fields. A proof that 
Eqs. (16) and (17) fulfill Eq. (5) can be found in Kossecka [44]. Here, Gij is the elastody- 
namic Green tensor of the anisotropic Navier equation defined by 

[kk f> dtt - Cijkldjdl] Gkrn = Sirn ^ {t)S{r) , (18) 

where 6{.) denotes the Dirac delta function and 6ij is the Kronecker delta. For isotropic 
materials, the tensor of elastic moduli reduces to 

Cijki = A 5ij5ki + /^{SikSji + SiiSjk) , (19) 

where A and are the Lame constants. Substituting Eq. (19) in Eq. (18), the isotropic 
Navier equation for the dynamic Green tensor was obtained 

[Sikpdtt - Siki^A- {X + 1^) didk] Gkm = SimS{t)5{r) , (20) 

where A denotes the Laplacian. When the material is isotropic and infinitely extended, 
the three-dimensional elastodynamic Green tensor reads [47, 48, 2, 3] 



+ rSit-r)drl (21) 



where r = x"^ + y"^ + z^. It should be pointed out that Eq. (21) is the retarded Green 
tensor. Here, cl and ct denote the velocities of the longitudinal and transversal elastic 
waves (sometimes called P- and S- waves) . The two sound velocities can be given in terms 
of the Lame constants (ct < cl) 



The elastodynamic Green tensor (21) is a tensor- valued distribution with support along 
the two sound cones r = Ct^ and r = c\^t as well as in the region between them. Eq. (21) 
consists of near-field and far-field terms. The first terms in (21) decay as 1/r and thus, 
they are the far- field terms. The last term in (21) decays more rapidly hke 1/r^ which 
gives the near-field term (see, e.g., [3]). 
If we use the relation 



2 rrjc-Y /"I/ct 

- / T(5(i-T)dT= / K?>{t-Kr)AK, (23) 
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substitute the Green tensor (21) into Eqs. (16) and (17) and perform the integration in 
t', the retarded elastic fields produced by an arbitrary three-dimensional distribution of 
dislocations are found 







Ri Rj \ 






i?3 ) 



Cl rC 



and 



-)^^^«/im(r',yd«;|dr', (25) 



where the so-called retarded times are given by 

R 

t^^t , (26) 

Ct 

iL = i--, (27) 
Cl 

t^^t-KR, (28) 

and K is a dummy variable with the dimension 1/ [velocity]. Here and are the 
transversal retarded time and the longitudinal retarded time, respectively. The retarded 
time is an effective retarded time for the integration with the limits (1/cl, 1/ ct) ■ Since 
Cl > Ct, the retarded times fulfill: tx > ^l and G [^l^^t]- Because the integrals (24) 
and (25) are evaluated at the retarded times, they are called retarded elastic fields. Here 
R = \r — r'\ is the distance from the source point r' to the field point r and is independent 
of t. The retarded elastic fields (24) and (25) are just integrals in r'. Here V is a volume 
integral in W. The sources T and / at the position r' depend on the retarded times. 
The retarded elastic fields at the position r and time t contain contributions from the 
past sound cones. Elastodynamic fields and waves propagate with finite velocities. Thus, 
there always is a time delay before a change in elastodynamic conditions initiated at 
a point of space can produce an effect at any other point of space. This time delay is 
called elastodynamic retardation. The retarded elastic fields (24) and (25) consist of three 
characteristic parts. The first term is the transversal one, transmitting with speed Ct, and 
it corresponds to S'-wave motion. The second term is the longitudinal one, transmitting 
with speed cl, and it corresponds to P-wave motion. Finally, the third term is neither 
longitudinal nor transversal and it gives contribution arriving at the speeds between the 
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two characteristic ones, which shows that this factor represents a combination of P-wave 
and 5'-wave motion. For a three-dimensional distribution of dislocations, we conclude that 
dislocations are retarded but not haunted by its past in contrast to a straight dislocation. 
Now we carry out the differentiations in Eqs. (24) and (25) and use the relations 

Rk R 

dkTlmn{r\ tret) = ^ dtTimn{r' , tret) , dkkmir' , ^ret) = 7r dtllm{r' , ^ret) , 

cK cK 

with tret = *T,*L,*«, C = Ct, Cl, 1/k , (29) 

in order to find the retarded elastic field more explicitly. We obtain 

l ^ f \ ^ I ^ijRk + SjkRi + ^ikRj 2)RiRjRk 



1 f ^ RjRi\ Rh „ „ , , , 1 / 5iuRj -\- SjifRi 3RjRiRh\ „ , , , 

- ^ ^*^^-(^ ' ^t) - ^ V,3 ' - J Timn{r\ t^) 



1 RiRjRk r^rr. /„/ + \ f SijRk + 3SjkRi + 3dikRj 9RiRjRk\ f '""^ („i + \a^ 

? P^ O't-'imnl^ , IlJ — ^ j J J-lmnK^ il'Kj'^K' 

+ (^-^ - ^ 9tr,^„(r', t,) d«| dr' 



+ ('^ - ^) ''a<i"j,n('-',(Jd« ^ dr' (30) 



and 



1 ^ /" f 1 / (^ij-Rfe + <^jfe-Ri + <^ife-Rj 2>RiRjRk 



-R^jR^ c^t/^.n(r- , tT) - ^ ^ J Iimir , ^l) 



1 RiRjRk T ^ / J. \ f SijRk + "^SikRi + iSikRj 9RiRjRk\ [^^'^'^ . / / . \ i 
+ — ^4 ^t^;m(r-', ^l) - ( ^ '-^^ ] / « /;m(r-', t^) 



These are the general retarded elastic fields produced by a time-dependent dislocation 
distribution. We now see from Eqs. (30) and (31) that the elastic distortion has three 
sources: the dislocation density T, the time derivative of T, and the time derivative of 
I. And we see that the elastic velocity has two sources: the dislocation current J, and 
the time derivative of I. All of these sources are retarded due to the retarded times. 
Equations (30) and (31) show that, since both equations contain the time derivative of /, 
the elastic distortion and the elastic velocity fields are created by the same time- variable 
dislocation current depending on the retarded times. Thus, the time-changing of the 
dislocation current / is the common source for /3 and v. The three (5ij-terms multiplied 
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by l/cx-factors in Eq. (30) and the two Sij-terms multiplied by l/cx-factors in Eq. (31) 
are analogous to the Jefimenko formulae for the electromagnetic field strengths (see [49, 
50, 23, 51]). Originally, Jefimenko [49] (see also [50]) derived the proper time-dependent 
generahzations of the Coulomb law and the Biot-Savart law as causal solutions of the 
Maxwell equations^. It is obvious that the retarded elastic fields (30) and (31) are more 
complicated than the Jefimenko equations in electrodynamics due to the tensor structure 
of the elastodynamical Green tensor and the dislocation-source fields. Nevertheless, we 
may call Eqs. (30) and (31) the Jefimenko type equations of dislocations. The causal 
dependencies of the elastodynamic phenomena of the motion of dislocations are described 
by Eqs. (30) and (31) which are exact solutions of the Navier equations (14) and (15) 
involving integrals over retarded dislocation sources. 

Performing the /t-integration and using Eqs. (19) and (22) the static limit of Eq. (30), 
which gives the (static) dislocation version of the Biot-Savart law is calculated as 

(32) 

It gives the correct expression for the elastic distortion tensor in terms of the dislocation 
density tensor T. Eq. (32) is the explicite expression of the isotropic version of the so- 
called Mura- Willis formula (see, e.g., [43, 52]). Here u is the Poisson ratio with A — 
2//i//(l - 2u) and u = A/[2(A + //)]. 



4 A non-uniformly moving dislocation loop 

Investigating the non-uniform motion of a dislocation loop, we consider a closed loop 
of arbitrary shape (planar or non-planar) that moves arbitrary. The dislocation density 
tensor and the dislocation current tensor of a dislocation loop at the position s{t) are 
represented by fine integrals of the form (e.g. [16, 44, 17]) 

Tijkir, t) = bi€jki (f S{r- s{t)) dLi{s{t)) , (33) 

JL{t) 

Iij{r,t)^kejki(f Vk{t)S{r-s{t))dLi{s{t)), (34) 

JL{t) 

where V = s denotes the velocity of the dislocation loop at any point s{t) on the loop, bi 
is the Burgers vector, L{t) is the dislocation loop curve at time t and dLi is a line element 



^In electrodynamics, the Jefimenko formulae for the electric field strength E and the magnetic field 
strength B are given by [50, 23]: 

47reo Jv \ R cR c R J 

-<-)^i^X(^^^^-^^^^^)--. 

where p is the electric charge density, J denotes the electric current density vector, c denotes the speed 
of light and eg is the permittivity of vacuum. 
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along the loop. Since L{t) may change its position with time t, it has a more complicated 
structure than in the static case. L(t) is the collection of all points on the dislocation line. 
Only subsonic source-speeds will be admitted (| V| < ct). If we substitute Eqs. (33) and 
(34) into (16) and (17) and after the integration in r', we obtain for the elastic fields of 
a dislocation loop 



Pim{r, t) = dk f emnp Cjkin Gij{r - s{t'),t - t') hi dLp{s{t')) dt' 

J -co JL{t') 

+ dt[ I pGij{r-s{t'),t-t')bjemnpVndLp{s{t'))dt' (35) 

J-oo JlW) 



and 



Vi{r,t)^dk f CjkimGij{r - s{t'),t-t')biemnpVndLp{s{t'))dt' , (36) 

J-oo JLit') 

which are in agreement with the formulae given by Markenscoff [19]. 

Integration in the expressions (35) and (36) may be performed in time. In this way, 
we find the elastic fields as hne integrals around the loop L{t') 

Pimir, t)^dk f emnp Cjkin (f>ij k dLp{s{t')) + dt f p Aijk bjemkp dLp{s{t')) (37) 

JL{t') JLit') 

and 

Vi{r, t)^dk f Cjkim Ajn biemnp dLp{s{t')) , (38) 

JLit') 

where we have introduced the elastodynamic Lienard-Wiechert tensor potentials (pij and 
Aijk in terms of the elastodynamic Green tensor Gij{r — s{t'),t — t'): 

/t i*oo 
/ G,,(r -r',t- t') 5{r' - s{t')) dr' dt' 
-oo J —oo 

= r Gijir-sit'),t-t')dt', (39) 

J —OO 

/t POO 
/ Gij{r -r',t- t') Vk{t') S{r' - s{t')) dr' dt' 
-oo J —oo 

= r Gij(r-s(t'),t-t')Vk(t')dt'. (40) 

J —oo 

More precisely, Eqs. (39) and (40) are the elastodynamical Lienard-Wiechert tensor po- 
tentials corresponding to delta-point sources acting on the position s{t') which have to 
be integrated over the loop hne element dLp{s{t')) in Eqs. (37) and (38). The property 
of the Green tensor (18) leads to the following wave equations for the elastodynamic 
Lienard-Wiechert tensor potentials (39) and (40) 

[Sik p dtt - Cijkidjdi] (t)km = 5{r- s{t)) , (41) 
[5ik P dtt - Cijkidjdi] Akmn = 5im K(^) 5{r - s{t)) , (42) 
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with Eq. (19). 

Substituting the elastodynamic Green tensor (21) into Eqs. (39) and (40), the inte- 
gration in time may be performed. Before the integration R = \r — r'\ is a function of r 
and r'; after the integration, which fixes r' — s{t), R — \r — s{t)\ is a function of r and 
s{t). The structure of the Green tensor (21) produces three characteristic integrals which 
we have to calculate. We express the integrals in terms of retarded variables by appeal 
to the relation [53, 54, 13] 



6{f{t'))g{t')dt' 



git') 



\df/dt'\ 



(43) 



at/(t')=0 



Mathematically, the factor l/|d//dt'| is the Jacobian of the transformation from t' to the 
new integration variable f{t'). This mapping between the two variables is one-to-one if 
the Jacobian is different from zero. A sufficient condition for this is that the velocity of 
the source (dislocation) is less than the wave speed (subsonic motion). The first integral 
can be carried out with 



5{t-t' -\r-s{t')\/cT) 
\r-s{t')\ 



dt' 



R-R-V/ct 



(44) 



t'=tT 



where R = r — s{t') and V = V{t'). R is the distance vector from the position of the 
source s, the sender of elastic waves, to the point of the observer r, the receiver of the 
elastic waves. The second integral is 



/ 



(r - sjt'Mr - sit% 
\r-s{t')\^ 



5{t - 1' 



\r-s{t')\/c^,^)dt'^^'^^ 



R-R-V/c^T 



(45) 



Using the relation (23), we perform the third integral as follows 



/( 



3(r - sjt'Mr - s{t')), 
\r-s{t')\^ 



\r - s{t')\ 

K dK 



1/cT 



n5{t-t' - n\r - s{t')\) dn dt' 



I/CL 



R-kR-V 



(46) 



From the argument that the delta functions vanish f{t') = in Eqs. (44)-(46), we obtain 
the condition 



t-t' - \r - s{t')\/c^O, 



with 



c = Ct,Cl, 1/k. 



(47) 



Unfortunately, the retarded times tc — t'{r, t) are not given directly, and they must be 
determined by solving Eq. (47) what can be quite tedious. Only in some simple cases 
tc is easy to find. If the dislocation loop is moving with subsonic speed, the solution of 
Eq. (47) is unique. The retarded times are a result of the finite speeds of propagation for 
elastodynamic waves. In Eqs. (44)-(46) we have used the relation 



d/(t') 



dt' 



= 1 



R V 



t'=tret 



cR 



t'=tret 



>0 for|y|<CT, c = ct,cl,1/k and t^et = ^t, ^l, , 

(48) 
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where f(f') = t-t' -\r- s{t')\/c. 

Carrying out the t'-integration in Eqs. (39) and (40), we find the exphcit expressions 
for the elastodynamic Lienard-Wiechert tensor potentials of a 'point dislocation source' 
acting on the position s{t') 



4:np(f)ij{r,t) = — 



+ 



I/CT 
1/cL 



i?2 JR-R-V/CT 



t'=tT 



+ — 

Ct. 



RiRj 1 
i?2 R-R-V/cl 



3RiRi 



i?2 



K, d/t 



R-kR-V 



4:npAijk{r,t) = — 



RjRi 



i?2 



V, 



+ 



i 



1/cT 



R- R-V/ct_ 







1 




t'=tT 





RiRn 



i?-i2.\//cL 



3 RiRj 



t'=tL 

(49) 

t'=tL 

(50) 



where R and V arc to be evaluated at the corresponding retarded times. The elastody- 
namic retarded potentials (49) and (50) fulfill 



Ajk{r, t) = Vfc(tret) 0y (y, t) . 



(51) 



The first 5jj-terms in Eqs. (49) and (50) have an analogous form as the original Lienard- 
Wiechert potentials of a point charge in the electromagnetic theory'' (see, e.g., [22, 23]). 
The Lienard-Wiechert tensor potential (49) is analogous to the displacement vector, which 
is the Lienard-Wiechert vector potential, of a non-uniformly moving point force in elas- 
todynamics found by Lazar [56]. Due to the appearance of two velocities of the elastic 
waves, the elastodynamic Lienard-Wiechert tensor potentials have a more complicated 
but rather straightforward structure. Like the retarded elastic fields (24) and (25), the 
Lienard-Wiechert tensor potentials (49) and (50) consist of three characteristic terms, 
again the first term corresponds to the S-wsve motion, the second term represents the P- 
wave motion, and the third term corresponds to a combination of P and S motion. Thus, 
at the retarded positions the dislocation loop is the source of S-, P- and the combination 
of P- and 5*- waves. 

Finally, if we substitute Eqs. (49) and (50) in Eqs. (37) and (38), the elastic fields read 



''The original Lienard-Wiechert potentials (scalar potential </> and vector potential A) of a point charge 
read [55, 50]: 



R{t')-R{t')-V{t')/c 



A{r,t) 



V{t') 



R{t')-R{t')-V{t')/c 



where q is the electric charge. Here tc denotes the retarded time with respect to the velocity of light. 
They fulfill the Lorentz gauge condition: -|- c^divA = 0. 
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in terms of the Lienard-Wiechert tensor potentials 



Pim{r,t) 



Cjkln ^mnp ) 2 



47rp 



1 

+ — 



RiRj 



L{t') 
1 



ut^. R^ R-R-V/cj^ 



JR-R-V/CT 

dLp{s{t')) 



t'=tT 



t'=tL 



1/cT 



+ 

'Lit') 



SRiRj 



K 



L{t') 



' R-kR-V 

RjRi \ Vk 



dK 



i?2 JR-R.V/CT 



t'=tT 



+ 



+ 



1 



RiRj 



R' R-R-V/cj^ 



I/CT 



1/cL 



L(t') 



3RiRj 



and 



Vi{r,t) = 



Cjklm bl (^mnp dk ^ ^ 



+ - 



RiRj 



L{t') 



dLp(s(t')) 

ft: Vfc 
R-kR-V 



RiRj 



dLp{s{t')) 



dn 



(52) 



'I] 



L(t') 



JR-R-V/ct 

dLp{sit')) 



dLpisit')) 



t'=t'T 



+ 



SRiRj 



R-kR-V 



dLp{s{t')) 



(53) 



where again R — r — s{t') and V = s{t') are to be evaluated at the corresponding 
retarded times which must be determined by solving Eq. (47). The fields (52) and (53) 
must be evaluated at some earlier times t' (the retarded times) and for the corresponding 
point s{t') on the dislocation loop. Also the line element dLp depends at every point 
s{t') on the retarded times. Thus, the elastic fields are given as line integrals around 
the dislocation line in terms of the Lienard-Wiechert tensor potentials of point sources. 
Due to the retarded times in Eqs. (52) and (53), L(t') depends on the retarded times and 
three curves L{tT^), L{ti^) and L^t^) appear in Eqs. (52) and (53). Because of the explicit 
dependence of the retarded times on s{t'), every point s{t') on the loop L{t') depends on 
its own retarded time. For a dislocation loop, the time dependence of the elastic fields is 
based on a retardation due to the retarded times which are functions of the variables r, 
s, and t. The elastic fields (52) and (53) at the point r and at time t receive contribution 
from every point s(tret) on the moving loop sending the elastic waves (S',P, mixed waves) 
at the retarded times (tx, ^l, tn)- For a moving loop these times are different not only due 
to different points on the loop, but also because the loop moves. The evaluation of the 
fields (52) and (53) is far from a trivial task. Thus, Eqs. (52) and (53) are complicated 
line integrals depending on the retarded times, and time and spatial derivatives outside 
the line integrals. 

Alternatively, the substitution of Eqs. (33) and (34) into the retarded elastic fields (24) 
and (25) gives directly the fields (52) and (53). Thus, the tensor structure of (52) and (53) 
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is inherited from the tensor structure of the retarded elastic fields (24) and (25). Because 
the time- variable of the sources T and I in Eqs. (24) and (25) is the retarded time tret, 
consequently the closed loop curve depends on the retarded times L(tret) in Eqs. (52) and 
(53). Although the structure of Eqs. (52) and (53) is elegant from the mathematical point 
of view, the evaluation of such expressions for a non-uniformly moving dislocation loop is 
very complicated due to the dependence of the retarded times. But this is the price we 
have to pay if we perform the integration in time and use the three-dimensional Green 
tensor (21) in order to obtain the elastodynamic Lienard-Wiechert tensor potentials. Since 
the relation of the retarded position to the present position of the dislocation loop is not 
in general, known, the elastic fields and the Lienard-Wiechert tensor potentials ordinarily 
permit only the evaluation of the fields in terms of the retarded positions and velocities 
of the dislocation loop. The complexity of the fields (52) and (53) is hidden behind the 
elegance of these formulae. 

On the other hand, if we substitute Eqs. (33) and (34) into the Jefimenko type for- 
mulae (30) and (31) and carry out the integration in r', we obtain 



1 r 1 



^ijRk + ^jkRi + ^ikRj 

SRiRjRk 



i?4 J R-R- V/ct 

RiRj \ Rk 1 



dLp(s(0) 



t'=tT 



dLp{s{t')) 



-1/cT 



J lie,. I Jut' 



V i?2 ; RR-R-V/CT 

I SjkRi + ^ikRj _ ^RiRjRk\ 1 

1^ R^ R^J R-R-V/ci^ 

RiRjRk 1 
^,,) R-R-V/ci^ 

^ijRk + ^SjkRi + 35ikRj QRiRjRk \ k, 



t'=tT 

dLp{s{t')) 



t'=ti. 



dLp{s{t')) 



I/CL lJL{t') 

+ dt / 



R-kR-V 



dLp{s{t')) 



Lit') 



SRiRj \Rk K 

)^ R- kR-V 



R^ 



1 



Lit') 

RiRj 



RiRj 

Sii 



Vk 



i?2 JR-R-V/CT 



dLp{s{t')) 
dLp{s{t')) 



t'=tT 



L{t') 



V) R-R-V/cl 



dLp{s{t')) 



t'=t-L 



3RiRj \ KVk , , ... 



dn 



(54) 
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and 

Vi{r,t) 



4:11 P 



mnp 



L{t') 



^ijRk + ^jkR-i + ^ikRj 



i?2 

3RiRjRk 



L(t') 



s,. 



RiRj \ Rk 



JR-R-V/CT 



dL,{s{t')) 



t'=tT 



i?2 J R R-R. V/ct 



dLpis{t')) 



Lit') 



SjkRi + ^ifc-Rj SRiRjRk 
Vr 



i?2 

RiRjRk 



Lit') 



jR-R-V/cj^ 
dLp(s{t')) 



t'=tT 

dLp{s{t')) 



t'=ti^ 



I/CT 
1/cL 

+ 9* 



L(t') 

1/cT 
1/cL 



(^ii-Rfc + "iSjkRi + 36ikRj 9RiRjRk 



L{t') 



i?2 

3 RiRj 



(5. 



i?4 



R-kR-V 

dLp{s{t')) 



dLpis{t')) 



t'=tn 



(55) 



Eqs. (54) and (55) are analogous to the Heaviside-Feynman formulae [57, 58] of elec- 
tromagnetic theory and we may call them the Heaviside-Feynman type formulae for a 
dislocation loop. They are more complicated than the original Heaviside-Feynman for- 
mulae [57, 58] (see also [59, 23, 51, 60]) for a point charge in electrodynamics^. Eqs. (54) 
and (55) contain terms without derivatives and terms with time-derivatives. They still 
possess a quite complicated tensor structure and line integrals depending on the retarded 
times. The causal dependencies of the elastic fields of a non-uniformly moving dislocation 
loop are described by retarded fine integrals along the loop. 



5 Static limit of a dislocation loop 

In this section, we give the static limit of the Lienard-Wiechert tensor potentials and of 
the elastic distortion of a dislocation loop as a check of the above results. 

In order to carry out the static hmit, we set V = and s(t') = r' and we substitute 
Eq. (22) and A = 2/iz//(l — 2z/) into Eq. (49). If we perform the integration in k and 
arrange in proper order the appearing terms, we find 



G 



(56) 



^In electrodynamics, the Heaviside-Feynman formulae for the electric field strength E and the magnetic 
field strength B of a non-uniformly moving point charge are of the form [51, 23]: 



E{r,t) = 



47reo 

q 



R 



E?{R-R-V/c) 
V xR 



c 



R 



R'^{R-R-V/c) 



1 



R{R-R-V/c) 
V xR 



R-R-V/c 



tc 



R{R-R-V/c) 
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with 

G,, = _ [2(1 - v)5,A - R , (57) 

where R = r — r' . Eq. (57) is the static three-dimensional Green tensor of the Navier 
equation [61, 1]. That means that the Lienard-Wiechert tensor potential reduces to 
the static elastic Green tensor in the static limit. In this manner, we find from Eq. (37) 
the so-called Mura formula (e.g. [15, 43, 62]) as static limit 



^ GmnpCjklnGij^kbl d-L'p , (58) 



which gives the elastic distortion of a static dislocation loop. 
On the other hand, the static limit of Eq. (54) is obtained as 

l^imir) = -^^^(Y^ Wp|(l - 2u) -3 + \ dL; , (59) 

which is the explicit form of Eq. (58). 



6 Straight dislocations 

We derive now the two-dimensional Lienard-Wiechert tensor potentials and the elastic 
fields of straight dislocations in our general framework. We choose the dislocation line 
iz parallel to the z-axis. The dislocation density and dislocation current tensors of edge 
dislocations are given by 

Tijk = bi ejkz iz5{R), lij = h ejkzVk ^z5{R). (60) 

For a screw dislocation, the dislocation density and dislocation current tensors read 

Tzjk = bz ejkz 4 5{R) , Izj = bz e^kzVk 4 > (61) 

where R = r — s{t) e and k — x,y. The tensors (60) and (61) fulfill, respectively, 
the relations 

hj — Tijk , Izj — Vk Tzjk ■ (62) 

Eq. (62) gives a relation between the dislocation current tensor, the dislocation density 
tensor and the dislocation velocity vector, which is valid for a single straight dislocation. 
Giinther [39] and Tcodosiu [36] derived such a relation between / and T for uniformly 
moving dislocations, that means for constant dislocation velocity V. If the relation (62) 
is valid, then / is a convection dislocation current. 

If the material is infinitely extended, the two-dimensional elastodynamic Green tensor 
of plane-strain reads [48, 63] 

(63) 
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and the elastodynamic Green tensor of anti-plane strain is given by 



1 H{t-r/cT) 



(64) 



where H{.) denotes the Heaviside step function and r = 

If we substitute Eqs. (63) and (60) in Eqs. (16) and (17) and perform the integration 
in r', we find for the elastic fields (no summation over z in Eqs. (65) and (66)) 



Pim{fi t) — ^mnz Cjkln dk(f>ij bi£z p dtAijk hjErnkz 



(65) 



and 



'^iiXi ^) — Gjklm 9k-^ijn ^l^mnz ^z 



(66) 



Using the definitions of the elastodynamic Lienard-Wiechert tensor potentials (39) and 
(40) and the Green tensor (63), we have here introduced the two-dimensional Lienard- 
Wiechert tensor potentials 



2np 



2np 



-00 
tT 

-ex; 



RiRj P 



+ 



RiRj 



Sij R^ 



RiRn 



R^ 



RiRj 



—00 

tT 



Vk{t') 
Vk{t') 



RiRj P 

RiRj ry 



+ 



R^ 



Sij R^ 



RiRi 



dt' 



+2 



Sij R^ 



RiRj 



% R' 



Si^jdt' 



(67) 



(68) 



The notation here is 

Ri — -^i ^i(^ ) ) 



t^t-t' , 



R' 



r,2 ' 



Si 



R^ 



(69) 



The two retarded times tx = t' and tL = t' are the roots of 5*1. = and S'l = 0, respectively, 
which are less than t. The solving of the conditions = and 5'l = is for a general 
motion non-trivial and can be very complicated. For subsonic dislocations the solutions 
for the retarded times tT and are unique. The general expressions (65)-(68) contain 
the fields for gliding as well as climbing dislocations. If V\\b, the expressions (65)-(68) 
describe a gliding edge dislocation and ifV -L b, Eqs. (65)-(68) give the fields of a climbing 
edge dislocation (see, e.g., [8]). For these cases, we can recover from Eqs. (65)-(68) the 
explicit expressions given by Kiusalaas and Mura [6], Lardner [1] and Lazar [8]. 

For a screw dislocation, if we insert Eqs. (64) and (61) in Eqs. (16) and (17), and 
performing the integration in r', we obtain for the elastic fields 



f^zmij' 1 ^) ^mnz ^zkzn 9k4'zz ^z ^z ~^ P ^t-^zzk ^z^mkz 



(70) 
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and 

Using the definitions Eqs. (49) and (50) and the Green function (64), the two-dimensional 
Lienard-Wiechert potentials of anti-plane strain are 

1 Z"*"^ 1 1 Vi.(f') 

<l^zz{r, t) = / — dt' , Azzk{r, t) = / ^ dt' . (72) 

It should be noted that the index z in Eqs. (70)-(72) is a fixed index (no summation 
over z) due to the reduction from 3D to the anti-plane strain problem. Prom the general 
expressions (70)-(72), we can reduce all the field components given by Eshelby [5], Nabarro 
[64], Lardner [1], and Lazar [8]. 

The two-dimensional Lienard-Wiechert potentials (67), (68), and (72) are time-integrals 
over the history of the motion and thus, they possess an afterglow. For that reason, a 
straight dislocation is haunted by its past as Eshelby [4] mentioned. Nevertheless, for the 
two-dimensional Lienard-Wiechert potentials (67), (68) and (72) we are left to evaluate 
time-integrals of considerable complexity which only in some simple cases yield results of 
elementary functions in a closed form. Also the calculation of the retarded times is not a 
trivial task. For straight dislocations, the static hmit of the elastic fields of non-uniformly 
moving straight dislocations was given in [8]. 

7 Near-field approximation of straight dislocations 

In the near-field approximation, it is important to determine the character of the sin- 
gularities of Eqs. (67), (68) and (72) at i? — > 0. Thus, we have to calculate only the 
near-field approximation of the integrals in Eqs. (67), (68) and (72). We follow the near- 
field approximation for cylindrical waves given by Courant and Hilbert [65], Whitham 
[66] and Barton [13]. We consider here dislocations which arc at rest and at t' = they 
start to move. This problem can be viewed as the superposition of the static equilibrium 
problem for t < with the dynamic problem for t > 0. Since the wave propagation na- 
ture of the solution is of interest, we consider the superposition-related problem obtained 
by subtracting the equilibrium solution from the complete solution. Therefore, the lower 
integration limits of the Lienard-Wiechert potentials (67), (68) and (72) are changed from 
-00 to (see also [67, 68, 69, 21]). 

Wc start with the near-field approximation for a screw dislocation. In the near-field 
approximation, the retardation is negligible. Thus, we have V{t') ~ V{t) and R ~ 
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r — s{t') ^ r — s{t). The two-dimensional Lienard-Wiechert potentials (72) reduce to 

2 , , H(t- R/ct) , 

27rpc| 0,,(r, t) ^ 1^ dt', r = t - t' 

dr 













- i?/cT) 







dr 

H{t - R/ct) In 



c^t / c^t^ 



« (« - B/ct) In ^ + O(R^) (73) 
and A22fe(r,i) — Vk{t)(f)zz{r,t). For i ^ R/ct, Eq. (73) simplifies to (see also [13]) 

If we substitute (74) into Eqs. (70) and (71), we obtain the singular terms of the elastic 
fields of a moving screw dislocation 

Pzm{r,t) = -€rnnz^;^^[C;,kzn^ - pVn{t) pVn{t)ln^-- j (75) 



and 



( Mir ^^^n{t) . . 

Exphcitly, they read 

R (77^ 

Pz.[r,t)- J^^ ^ ^2 4 it: j 

R (.f\ ^^^^//^l ^'^^)\^^ V^myjt) Ry V;(Oi 2cTn 

«r,t) = ^|(^i-^j^ — 4^:r^"^^"^/ ^^^^ 

..(.,t)^^W)_^^. (79) 

It is obvious in Eqs. (77)-(79) that a screw dislocation contains a singularity and 
a logarithmic singularity in the near-field. The acceleration terms give the logarithmic 
singularity. 

Also Eqs. (77) and (78) contain the correct static limit with V = and s{t) — constant 
for the clastic distortion of a screw dislocation with the Burgers vector b — (0, 0, hz) given 
by deWit [70] 
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For convenience we have chosen s = 0. 

For the near-fields of edge dislocations, we need the following integral approximations 



H{t - R/c) 



r h{t - . 

Jo 



dt' 



H{t - R/c) 



-Jo 




R^/c^ 


= H{t 


-R/c) 


f- 

J R/c y 


= H{t 


-R/c) 


(R^ 


^ H{t 


-R/c) 


\2c2 



dr 



In 
In 



ct 



cH"^ 



+ 



t , i?2 



0{R' 



(82) 



and 

ft 



[ H(r - R/c) - i?7c2 dt' = / H(r-R/c) ^r'^-R^/c^ dr 
Jo Jo 

^H{t-R/c) [ y/r^-Ryc'^dr 

J R/c 



-H{t-R/c){ - — In 



R^ 

i?2 



ct 

R^\W 



t I ^ i?2 



2ct e i?2 



Hit-R/c)[-^-X.^-^^--^-^Om). 



(83) 



Using Eqs. (82) and (83), we obtain for the Lienard-Wiechert tensor potentials (67) and 
(68) 



RR / R^ 
27rp<P^j{r,t) ^ H{t - R/c^){ -^tJt^ - 



Cht 

R 



cjf 



2i?2 



H{t-R/c,){^Jt^-^ 



C^t 



+ 



+ 



2i?2 



(84) 



and Aijk{r, t) — Vk{t) (j)ij{r, t). If t ^ R/ct and t » -R/cl, then (84) simplifies to 



<l>ij{r,t) 



1 \ RjR 



„2 2 



i?2 



2cTt 



2CLt 

R 



(85) 
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and Aijk{r,t) — Vk{t) (f)ij{r,t). We calculate 

1 



dk(l>ij{r,t) 



87r//(l - u) 

dkAijk{r,t) = Vk{t)dk(j)ij{r,t) and 



(86) 



dtAijk{r,t) = 



87r/x(l-i/) i?2 



■3 ^ yk{t) 



1 , 2cTt 1 , 2cLt 
In — h In 



R 



R 



+ 



87r/x(l - 



(3 - 4z/) 5, 



RV 



R,Vj+RjV^^^RiRjR-V 



^7) 



The first term in Eq. (87) is non-singular and, thus, it does not contribute to the singular 
near-field. Substituting (86) and (87) into Eqs. (65) and (66), we find for the singular 
near-fields of the elastic fields of edge dislocations 



Pim{r,t) 



87r/x(l - v) 



jkln 



pVn{t) 



(3 - 4z/) 5u 



RV 



Rk 

RiVi + RiVi 



(3 - 4Z/) 5ij -J^-"ik^- "jkj^2 



Oik- 



Ri 



i?2 



+ 2 



^jk^ + 2 

RjRiR-V 



RiRjRk 
R^ 



+ 



An 



1 , 2cTt , 1 , 2cLt 
In ^ m 



R 



R 



and 

Vi(r,t) = 



l^zCjklmyn{t) 

87r^(l - u) 



Rk 



Ri 



Rj 



/o A \ z "^k r- s I o Ri^j^k 

(3-4.)5,-^-5..^-<J,.^ + 2^^ 



(88) 



(89) 



where R = r — s{t). It can be seen in Eqs. (88) and (89) that edge dislocations contain 
a 1/R singularity and a logarithmic singularity in the near-field approximation. The 
acceleration terms in Eq. (88) give the logarithmic singularity. Eqs. (88) and (89) are 
valid for climbing and gliding edge dislocations. 

It is also important to note that Eq. (88) gives the correct static limit with V = 
and s{t) = constant for the elastic distortion of an edge dislocation with Burgers vector 
6 = (6^, 0, 0) given by deWit [70] 



47r(l - ly) 



(1 - 2i/) + 2 



X 



X 



47r(l — u) 



X 



47r(l - u) 

bx y_ 

47r(l - v) 



;i - 21/) + 2 ^ 

(3 - 2z/) - 2 - 



X 



2^) -2-. 



(90) 
(91) 
(92) 
(93) 



Again we have chosen s = 0. 
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It is important to mention that the history of the motion is not visible in the singu- 
larities of the non-uniformly moving dislocations (75), (76), (88), and (89). The singular 
terms are determined by the current situation. Both screw and edge dislocations are char- 
acterized by two types of singularities, namely, 1/R and logarithmic singularities. The 
logarithmic singularity is characteristic for the acceleration motion. The logarithmic sin- 
gularities (acceleration terms) appear only in the elastic distortions (75) and (88) and not 
in the elastic velocity terms (76) and (89). The character of the singularities of the near 
fields is in agreement with the results given by Markenscoff [71], Callias et al. [72] and Ni 
and Markenscoff [73, 74]. A near- field approximation of singular integrals giving simple 
coefficients of the singular terms has been used. In this simple near-field approximation, 
we have neglected the retardation. 

8 The Mura dislocation tensor potentials, the re- 
tarded dislocation tensor potentials and the Lienard- 
Wiechert tensor potentials of a dislocation loop 

On the other hand, Mura [75, 15] has written the elastic fields in terms of so-called 
dislocation tensor potentials in the following way 

Vi{y-,i^ — ^mnpC jklm 9k-^ijnlp ■ (^^) 

The dislocation tensor potentials (f)ijki and Aijkmn, introduced by Mura [75, 15], are defined 
by (see also [76]) 

/t poo 
/ Gij{r-r',t-t')aki{r',t')cir'cit' , (96) 
OO J — OO 

/t poo 
/ Gij{r -r',t- t') Vkmnir', t') dr' dt' , (97) 
-oo J —oo 

and they fulfill the following wave equations 

[^ik P dtt — Cijkldjdlj (pkmpq = ^im Oipq , (98) 
■^kmnpq ^im ^npq i (99) 

where Vkmn is Mura's dislocation velocity tensor [14, 77]. Although Vkmn possesses 27 
components, only 9 components which are skew-symmetric in some indices {Vkmn = —Vnmk 
and Aijkmn = —Aij^mk) enter Eqs. (94) and (95). For a straight dislocation with Burgers 
vector bm, dislocation line direction i„ and dislocation velocity Vk, it reads: Vkmn ~ 
Vkbm'^n^{R) — VkOimn- Here CKy denotes the usual dislocation density tensor (10). The 
dislocation current tensor (7) can be written in terms of Mura's dislocation velocity tensor 
as follows 



hj = ^jkiVkii ■ (100) 
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Substituting Eqs. (10) and (100) into Eq. (9), the Bianchi identity (9) reads 

oiii + Vkii,k - Viik,k = , (101) 

which leads to the foUowing 'gauge' or side condition for the dislocation tensor potentials 
(see also [15]) 



Hjkl 



~l~ -^ijmkljm ■^ijlkm,m . (102) 



If we substitute the Green tensor (21) into Eqs. (96) and (97) and perform the inte- 
gration in time t\ we find the retarded dislocation tensor potentials^ 



and 



Hjkmn 



+ - ^) t Mr', Q d«:| dr' (103) 

+ - ^) j^''\v,„(r',t,)dKjdr', (104) 



where the retarded times ^t, and are given by Eqs. (26)-(28). Because the integrands 
are evaluated at the retarded times, these fields are called retarded dislocation tensor 
potentials. They arc the causal sohitions of the inhomogencous Navicr equations (98) 
and (99). The retarded dislocation tensor potentials (103) and (104) satisfy the gauge 
condition (102), this can be easily checked. In the static case, the retarded dislocation 
tensor potentials (103) and (104) reduce to 



'^-'='(^)^16vrMl-.)/v^ 



(3-4.)5,. + ^ 



aki{r')dr', (105) 



and Aijkmn{f) = 0. On the other hand, if we substitute Eqs. (103) and (104) into 
(94) and (95) and use the relations (100), (11), and (25), we obtain the Jefimenko type 
equations (30) and (31). 

If we compare Eqs. (94) and (95) with (37) and (38), we find 

(l)ijmn{r,t) = f (l)ijbmdLn{s{t')) , Aijkmn{r,t) ^ (t Ajfe ^'m dL„(s(i')) , (106) 

JL{t') JL{t') 



®In electrodynamics, the retarded electromagnetic potentials were originally introduced by Lorenz [78] 
and they read [50, 23]: 

V, 1 f Pir',t-R/c) 1 f J{r',t-R/c) 

where p is the electric charge density and J denotes the electric current density vector. The idea of a 
retarded scalar potential was first developed by Lorenz [79] in 1861 when studying waves in the theory of 
elasticity. The retarded potentials fulfill the Lorentz gauge condition: + c^divA = (see, e.g., [50, 49]). 
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that means that the dislocation tensor potentials are given as line integrals of the three- 
dimensional elastodynamic Lienard-Wiechert tensor potentials of a so-called 'point dis- 
location source' (49) and (50) along the dislocation line L{t'). The dislocation tensor 
potentials of a dislocation loop (106) may be identified with the actual Lienard-Wiechert 
tensor potentials of a dislocation loop. If we substitute Eqs. (49) and (50) into (106), they 
read explicitly 
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47rp \ ( 
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(107) 



and 
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(108) 

where V = V{t'), R = R{t'), s{t') and L{t') have to be computed at the corresponding 
retarded times. Thus, the Lienard-Wiechert tensor potentials of a dislocation loop (107) 
and (108) are line integrals of the Lienard-Wiechert tensor potentials of a dislocation 
point source (49) and (50) integrated over three loop curves L(tT), L{ti^), and Llt^) 
at the corresponding retarded times and with the integrands dL„(s(tT)), dLn{s{tj^)), 
and dL„(s(tK;)). Eqs. (107) and (108) are retarded line integrals. Thus, the Lienard- 
Wiechert tensor potentials of a dislocation loop are the line integrals of a 'dislocation 
point source' acting on s{t') integrated over the dislocation loop L{t') and evaluated at 
the corresponding retarded times. The elastic fields of a non-uniformly moving dislocation 
loop are obtained by substituting Eqs. (107) and (108) into (94) and (95) reproducing 
the formulae (52) and (53). The fields produced by a dislocation loop can therefore be 
computed by first determining the Lienard-Wiechert tensor potentials (107) and (108) and 
then obtaining the elastic fields by means of the relations (94) and (95). More directly, the 
Lienard-Wiechert tensor potentials of a dislocation loop (107) and (108) can be obtained 
by substituting the dislocation density tensor and Mura's dislocation velocity tensor of a 
dislocation loop 



o(ij{r,t) ^bi i S{r - s{t)) dLj{s{t)) , 

JLit) 

Vkij{r, t) = bi(f Vk{t) S{r - s{t)) dLj{s{t)) , 

JLit) 



(109) 
(110) 
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into Eqs. (96) and (97) or in Eqs. (103) and (104). Eqs. (109) and (110) can be read-off 
from Eqs. (33) and (34), using the relations (11) and (100). 

In tlie static case, tlie Lienard-Wiechert tensor potentials of a dislocation loop (107) 
and (108) reduce to 



(3 - 4.)<^., + ^ 



di^;, (111) 



and Ajjfem„(r) = 0, where R = r — r'. 

Comparing Eqs. (94) and (95) with (65) and (66), we obtain for the two-dimensional 
Lienard-Wiechert tensor potentials of plane strain 

4'ijmz 4'ijbm^z i -^ijkmz -^ijk^m^z i (H^) 

and with (70) and (71), we obtain for the two-dimensional Lienard-Wiechert tensor po- 
tentials of anti-plane strain 

(t'zzzz 4'zzbz^z ) A^zkzz Azzkt^z^z ■ (H"^) 

From the condition (102), we get 

4- + ^ijk,k = (114) 

and 

<i)zz + Azk,k^o. (115) 

A 'gauge' condition like (115) was already mentioned by Lardner [1]. 

Therefore, it is obvious that the elastodynamical Lienard-Wiechert tensor potentials 
can be obtained from Mura's dislocation tensor potentials. Especially, for a dislocation 

loop the Mura dislocation tensor potentials give directly the Lienard-Wiechert tensor 
potentials as retarded line integrals along loop curves computed at the corresponding 
retarded times. 

9 Conclusions 

In this work, we have investigated the fundamentals of the non-uniform motion of dislo- 
cations based on electromagnetic analogies. We have examined and solved to following 
items: 

• retarded elastic fields/ Jefimenko type equations for dislocation fields 

• retarded dislocation tensor potentials 

• 3-D Lienard-Wiechert tensor potentials of point dislocation sources 

• 3-D Lienard-Wiechert tensor potentials of a dislocation loop 

• Heaviside-Feynman type equations for a dislocation loop 
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• 2-D Lienard-Wiechert tensor potentials of straight dislocations 

• singularities of the near-fields of straight dislocations. 

This analysis provides the treatment of general dislocation motion (with inertia effects) 
in terms of retarded fields. We have introduced the elastodynamic Lienard-Wiechert 
tensor potentials as fundamental quantities for the non-uniform motion of dislocations. 
For the general motion of a dislocation loop, the solution is obtained in terms of Lienard- 
Wiechert tensor potentials depending on the retarded times and in a more closed form 
as line integrals (Heaviside-Feynman type formulae). The Jefimenko type and Heaviside- 
Feynman type equations contain only time derivatives and no spatial derivatives. Thus, 
for the general motion of a loop the solution is evaluated in a closed form. It is concluded 
that causal dependencies of the elastic fields and Lienard-Wiechert tensor potentials of a 
non-uniformly moving dislocation loop are described by retarded line integrals along the 
dislocation loop. The non-uniform motion of dislocations leads to retardation effects which 
are important in the interaction between dislocations. Thus, the interaction between 
dislocations during the dislocation motion is not instantaneous, without being limited by 
the speed of the sound waves. We think that these results will open a new window for 
a better understanding of the non-uniform motion of dislocation fields. Especially, they 
can be used in simulations of discrete dislocation dynamics. Retardation effects should 
be considered in dislocation dynamics simulations at high strain rates and they become 
important at extremely high frequencies. 
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